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Abstract. The properties of the fluorescence Fe K a emission lines of a sample of 38 quasars (QSOs) observed with 
XMM-Newton are studied. These objects are included in the optically selected sample from the Palomar-Green 
(PG) Bright Quasar Survey with an X-ray luminosity 1.3 x 10 43 <L2-iofc e v < 5.1 x 10 45 erg s _1 and z<1.72. For 
each object in the sample, we investigated the presence of both narrow and broad iron lines in detail. A total of 20 
out of the 38 QSOs show evidence of an Fe K a emission line with a narrow profile. The majority of the lines are 
consistent with an origin in low ionization material, which is likely to be located in the outer parts of the accretion 
disk, the molecular torus, and/or the Broad Line Region. The average properties of the narrow Fe K a emission 
line observed in the sample are similar to those of Seyfert type galaxies as inferred from recent XMM-Newton and 
Chandra studies. A broad line has been significantly detected in only three objects. Furthermore, we studied the 
relationship between the equivalent width (EW) of the iron line and the hard band X-ray luminosity for radio 
quiet quasars. The analysis indicates that no clear correlation between the strength of the line and the hard X-ray 
luminosity is present, and our results do not show compelling evidence for an anticorrelation between these two 
quantities, i.e. the so-called X-ray Baldwin effect. 
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1. Introduction 

The iron K Q emission line was one of the first atomic 
features discovered in X-ray spectra of active galaxies 
ijMushotzkv et al1ll978|h Its origin was promptly inter- 
preted as ph otoelectric absorption follow ed by fluorescence 
emission (see lRevnolds fe Nowakl (|Fo03) for a review) . The 
combination of fluorescent yield (oc Z 4 ) and the (rela- 
tively) high abundance of iron justifies the strength of such 
a feature with respect to lines emitted by other atomic 
species. The K a line actually consists of two different lines 
that correspond to the 2p— >ls transition and that, in the 
case of 'neutral' iron (i.e. Fel-XVII), have an energy of 
6.404 and 6.391 keV, respectively. 

However, the resolving power of earlier X-ray spec- 
trometers does no t allow for separ ating both fea- 
tures (but see also iKaspi et alJ |2002) for results with 
Chandra HETGS) , so they are usually observed (and con- 
sidered) as a single emission feature. The maximum en- 
ergy for this line is 6.97 keV corresponding to emission 
from H-like iron. 

The X-ray spect ral survey of Seyfert g alaxies car- 
ried out with Ginga l|Nandra fe Pounddll99l found this 
feature in the vast majority of these sources. In the 



same years, several theoretical works suggested an ori- 
gin for this feature in circumnuclear optically-thick ma- 
terial (i.e. the disk and/or the m olecular torus), together 



terial (i.e. the disk and/or the molecular torus), together 
with X-ray Compton reflection llGeorge fe Fabianl 1991 : 
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Matt. Perola. fe PirJl99lllKrolik. Madau. fe Zvckilll99 
Ghisellini. Haardt. fe Mattlll994|) 

The study of the iron K Q emission line became one of 
priorities of X-ray astronomy after the discovery of a line 
with a broadened and skew ed profile in the AS CA spec- 
trum of MCG-6-30-15 and l|Tanaka et al.lll995|) . Such a 
broad profile was interpreted as the combination of trans- 
verse Doppler shift, relativistic beaming, and gravita- 
tional redshift. It was the first direct evidence of emis- 
sion from an accretion disk extending down to very few 
Schwarzschild radii from a supermassive blackhole. Iron 
lines are therefore a unique tool for studying strong grav- 
ity effects and mapping the innermost regions of the ac- 
cretion flow. 

The advent of XMM-Newton and Chandra has 
allowed more detailed investigations of this spectral 
feature and provided i mportant insights into it s prop - 
ertie s. In particular, lYaaoob fe Padmanabhar] l|2004h 
and IPounds fe Reeves! l)2002j) pointed out the ubiquity 
of an unresolved Fe K Q line in the spectra of bright 
Seyfert galaxies. Such a line should therefore emerge 
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from the outermost regions of the accretion disk, but 
non-disk contributi ons (i.e. from the Broad Line Region 
lYaaoob et all 120011) are also plausible. Co ntrary to pre- 
vious low-resolution ASCA observations ijNandra et alJ 
1997) iron emission lines with a broad profile appear to 
be quite rare in Chandra and XMM-Newton observations. 

It is worth noting that the Fe K a line has not 
commonly been observed in the X-ray spectra of high 
luminosity AGNs (i.e. QSOs) and, wh en detected, it is 
usually weaker than in Seyfert gala xies llReeves fe Turner 

mTo explain this effect, llwasawa fe Taniguchi 
suggested, based on Ginga data, the existence 
of an X-ray "Baldwin effect" such that, as the X-ray 
luminosity increases, the upper layers of the disk become 
progressively more ionized until any reproc essed feature 
is suppressed llRoss. Fabian, fe Yound ll999'). The X-ray 
Baldwin effect was also observed with ASCA (Nandra et 
al. 1997) and more recently with XMM-Newton (Page et 
al. 2004a). However, due to selection effects (i.e. relative 
faintness, lower density surface than in Seyfert galaxies) 
and to the limited spectral capabilities of past X-ray 
telescopes, knowledge about the properties of the iron 
emission line in QSOs is quite sparse. 

Our ext e nsive XMM-Newton program 

jPiconcelli et all l2004b|) is designed to fill this lack 
of information by detailed and uniform analysis of 40 
EPIC spectra of opticall y-bright QSOs selected f rom the 
Palomar-Green survey I Schmidt fc Greer] Il983h . With 
such a large sample we aim to provide useful insights 
into the characteristics of the Fe K a line and its origin. 
Furthermore, the homogeneity of our sample strongly 
prevents possible biases due to mixing different classes of 
AGNs (see Sect.©. 



2. XMM-Newton observation 

This study presents an analysis of the iron K band of 
all public observations of PG quasars available in the 
XMM-Newton Science Archive (XSA) as of February 
2004. Out of the 42 original objects, two (0003+199 and 
1426+015) were excluded due to large pile-up and two 
others (1001+054 and 1404+226) due to low statistics in 
the hard X-ray band, i.e. >2 keV. The spectra were ob- 
tained by processing the XMM-Newton observations us- 
ing the Science Analysis System (SAS) v. 5. 4.1. Sample 
properties and d ata red uction are described in detail in 
IPiconcelli et"ai] ll2004blh MOS1 and MOS2 0.6-10 keV 
spectra were combined and analyzed simultaneously with 
0.3-12 keV pn spectrum. 



3. Analysis of spectra 

3.1. Best fit models 

The spectrum of each QSO was analyzed using; 

XSPEC v.11.2, (jArna.udl Il99rfl . From IPiconcelli et"at) 



(2004a). we obtained the best fit model for the continuum 
emission of each QSO studied . 

Table \I\ shows the best fit model obtained for each 
source. A complete descripti on of these mode l s and t heir 
parameters are presented in IPiconcelli et alJ (j2004b) . In 
the following, all the energies of the QSOs are referred to 
the rest frame, and errors are given to the 90 per cent 
confidence level (i.e. A% 2 = 2.71) unless otherwise stated. 
Luminosities were calculated considering a flat cosmology 
with (Om^a) = (0 .3,0.7) and H = 70 km s" 1 Mpc -1 
l)Bennett et al.ll2003|h 

3.2. Spectral models for the Fe line 

First, we investigate the presence of fluorescent Fe emis- 
sion lines by adding an unresolv ed narrow Gauss i an line 1 
to the best fit models reported in IPiconcelli et~aT1 <|2004bj) . 
Whenever the addition of this component was not required 
by the data (i.e. its significance is < 95% according to an 
F— test), we then fixed the centroid of the line to 6.4 keV. 
In such way an estimate of upper limit to the rest-frame 
equivalent width (EW) of the neutral iron fluorescence 
emission was inferred. The best-fit values for the emis- 
sion line parameters are listed in Table Q An iron fluo- 
rescent emission line with a narrow profile is significantly 
detected, i.e. with an F — test criterion at, a significance 
level of >95%, in 20 QSOs ( 17 radio quiet and 3 radio 
loud). 

Fig. fjji shows the centroid of the iron fluorescence 
narrow emission line plotted as a function of 2-10 keV 
luminosity for the 20 QSOs with significant detection 
of such an emission feature (see Table QJ. There is no 
apparent trend with luminosity or radio-loudness. The 
centroid energy of the Fe K a line in all but four quasars 
(1115+407 1116+215 1244+026 1402 +261) correspond s 
to low ionization states, i.e. Fe I-XVII l|Makishimalll98fih . 
Fig. ^> shows the EW of Fe K-shell narrow emission lines 
as a function of luminosity in the hard band; however; 
however no obvious correlation emerges from this plot 
(see Sect. EH2~T|) . 

Besides the ubiquitous presence of a narrow Fe 
emission line centered at ~ 6.4 keV, recent works 
based on Seyfert 1 samples observed w ith both the 
Chandra llYaaoob fe Padmanabhanl l2004h and XMM- 
Newton l|Pounds fc Reeves! 120021) high signal-to-noise 
spectroscopic data have found in some objects evidence of 
additional underlying broad and skewed emission support- 
ing the original A SCA claims for a re lativistically broad- 
ened component (Na ndra et al.lll997lK Consequently, we 
added a disk line to the fitting model to account for 
any such "broad" component. For this we used the LAOR 
model following recent XMM-Newton results that ad- 
vocate a spinning black hole in nearby Seyfert galaxies 



The intrinsic width of the Gaussian line was forced 
to be the instrumental resolution at the given energy, 
FWHM~ 7000 km s" 1 for both pn and MOS, cr=0 keV. 



XMM PG Team: XMM-Newton observation of PG quasars 



3 



Table 1. Iron line properties. The table includes the name of the QSO, the best fit model of the data, the rest frame 
energy and intensity of the Fe K Q line modeled by a unresolved Gaussian line, the rest-frame equivalent width, the 
F-test probability, and the rest-frame equivalent width of a broad iron line modeled by a LAOR line. 



PG Name Model" L 2 -iofc e y 






Narrow Line 


Broad Line 






EW 


Intensity 


P(F - teat) EW 


(10 44 erg s" 1 ) 


(keV) 


(eV) 


(ph cm _2 s _1 ) 


(%) (eV) 



0007+106* 


A 


1.4 


6.4ft 


<60 


<6 x 10" 6 


<75 


220+"° (tt) 

w — 100 


0050+124 


E 1 


0.78 


6.47+°-?? 

— 0.11 


40 ± 30 


(4 + 2) x 10" 6 


97.5 


100 +4 ° (tt) 


0157+001 


A 


0.71 


6.4f 


< 330 


< 4 x 10 -6 


<75 


<1300 


0804+761 


D 


2.86 


6.38+0.06 


100+™ 


l-2ln ' 7 xl0 -5 


95.7 


250+?2° 

^ uyJ — 180 


0844+349 


E 


0.55 


6.4f 


< 100 


< 6 x 10" 6 


<75 


<1200 


0947+396 





2.26 


6 40+S'ii 

"■^ u -0.07 


120 ± 60 


3.1 ± 0.2 x 10~ 6 


99.8 


< 360 


0953+414 


E* 


5.41 


6.4f 


< 50 


< 4x 10~ 6 


<75 


<77 


1048+342 


E 


1.10 


6 34+° ^ 

u '°^— 0.04 


100+60 


1.9 + i'2 xl0~ 6 


97.6 


<160 


1100+772* 


E 


11.19 


6 40 +0 A° 

u,w -0.04 


60+40 


3+J x 10 -6 


95 


<110 


1114+445 


G 


1.45 


6 45+°-? 2 


100+5° 


3.3l 12 x 10~ 6 


99.99 


<33 


1115+080 


PL* 


65.3 


6.4f 


< 130 


< 2 x 10~ 6 


< 85 


1400+ 1 6 °° 


1115+407 


E 


0.85 


6.83+0.09 


130 + SS << 100 § ~) 


1.7to?x 10~ 6 


95 


<320 




C 


'i 1 


7 i y+0.07 


900+ 110 (<r 80 s ~1 


7+ 3 x 1 n~ 6 

— 5 


Qfi ^ 


500+ 1100 (tt) 


1 909-1-981 

_1_ i_ U |^ ij O -L 


c 


2.68 


6.4f 


<80 


< 2 x 10~ 6 


^75 

— 1 u 


<420 


1 9064-459 


PL 


14.9 


6.4f 


<350 


< 6 x 10 -6 


^75 


<580 


12114-143 


E 


0.50 


6 45 +0 - 2 


4 q+20 


17+1 1 x 10~ 6 


95.1 


OQQ + 300 


1 91 64-069 




4.8 


6.4f 


<r 70 


<1.6x 10~ 6 


^75 


<220 


12964-023* 


E 


51.05 


6.4f 


<7 


<6xl0 -6 


<75 


<17 


1 944+098 

± £'-±.'-±r\\j ZiU 




0.14 


6 65+ 07 
o.oo_ 18 


*3on 4- 9nn ( <r 1 ^n^ "i 


7+ 5 y in -6 






1 3074-085 

LOU i \\jO'J 


A 


1.19 


6.4f 


<110 


<3 xlO -6 


^75 


<920 


1 3094-355* 

-LOU£7 1 OOO 


A 


0.67 


a 07+0. 05 
u -°' -0.06 


1 80+ 80 


1 8 +0 ' 8 x 10" 6 


99.7 


<410 




p 

V . 


1 08 


fi 48+D HQ 


1 sn+i 1 n 


9 8+ 1 ' 6 v in -6 

^.0 X7 






lOJi | -LOO 




1.36 


6 44+0.04 

"■^-0.07 


1 50 + 80 

-LOU —1— OU 


4 ± 2 x 10~ 6 


97.7 


47n +i70 




( 

V 


1 A 
1 .-t 


7 47+O.O7 
'■ 4 '-0.3 


oqo -i- 190 ( mn§^ 

ZoU it 1ZU 1UU ) 


4+ 2 v in- 6 

^_3 A _1U 




<^ 1ZUU 


1407+265 


A 


41.20 


6.4f 


<30 


<1.0x 10 -6 


<75 


<66 


1411+442 


H 


0.25 


6 40+0.05 

°-^°-0.13 


2501"° 


2.9 + o^ x 10~ 6 


99.8 


900+300 


1415+451 


E 


0.4 


6.35 + 0.07 


110 ± 80 


1.5+1'? x 10" 6 


97.6 


7 qq+300 


1427+480 


C 


1.6 


6A0±° f 


qo+50 


1.3±g;l x lO" 6 


97.1 


<480 


1440+356 


D 


0.58 


6.4f 


< 80 


< 3 x 10~ 6 


<60 


<310 


1444+407 


D 


1.3 


6.4f 


<180 


< 2 x 10~ 6 


<75 


1QQQ + 5O0 

iuuu_ 900 


1501+106 


E 


0.49 


6.4f 


<60 


5< 1 x 10~ 5 


93 


<52 


1512+370* 


D 


8.84 


°- dz: -0.08 


120+60 


3.9I2 8 x 10~ 5 


99.5 


<240 


1613+658 


A 


1.78 


6.4f 


<60 


<4x 10~ 6 


<75 


<210 


1626+554 


D 


1.46 


6.4f 


< 160 


< 7x 10" 6 


<75 


<190 


1630+377 


D 


20.9 


6.48«;« 


600+300 


5+j x KT 6 


99.8 


<1400 


1634+706 


C 


127.9 


6.4f 


<82 


< 5 x 10" 6 


<75 


<680 


2214+139 


I 


0.48 


6.34+0.06 


80t 4 °, 


4.0+1.7 x 10~ 6 


99.9 


<150 


2303+029 


c 


16.2 


6.4f 


< 140 


< 1.6 x 10~ 6 


<20 


<160 



Radio loud objects. lt) In this source only an iron line with a broad profile was significantly detected. The detection of the 
LAOR line at 6.4 keV is significant, i.e. F-test > 95 %. ^ Best fit model: A: blackbody; B: multicolor blackbody; C: bremsstralung 
emission; D: power law; E: double blackbody; F: absori model ; G: double absori model + blackbody; H: partial-covering 
(cold) + Raymond-Smith; I: partial-covering(warm)+emission line; J: absorption + blackbody; PL: simple power law. Models 
with ^' include an additional cold absorption compone nt. All but one t ype of models (i.e. model PL) also include a power 
law component accounting for the hard X-ray band fsee IPiconcelli et alJ (|2004bh ). The EW labeled with correspond to the 
neutral iron line added in addition to the ionized one. 



Alms et al]l20niUFabiarJl2002l) . In our fit we assumed 
a fixed rest frame energy of the line to 6.4 keV, a power 
index of the radial emissivity law (3 equal to -2, an in- 
ner (outer) radius of 1.26(400) R g , and an inclination of 



30° 



The resulting EW of the broad component is listed 



2 Th ese are the same parameters used in lPorquet fe Reevesl 
(l2003h for the broad Fe K Q line observed in the XMM- 
Newton spectrum of the luminous quasar Q0056— 363 
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Table 2. Broad iron line properties. The table includes the name of the QSO for which a broad iron emission line is 
significant. The line was modeled by a LAOR line assuming a fixed rest frame energy of the line to 6.4 keV, an inner 
(outer) radius of 1.26(400) R g , an inclination of 30°, and (3 equal to -2 (see text for further details). The intensity and 
rest-frame equivalent width of the broad line, as well as the F-test probability, are also given. 



PG Name 


Intensity 
ph cm _2 s -1 


eV 


P (F-test) 

% 


0007+106 
0050+124 
1116+215 


(1.5 + 0.7) x 10 5 
(6 + 3) x 10" 6 
l.7t\l x 10" 5 


100t« 
500™ 


99.2 
99.6 
97.6 



I 



Lum 2-10 keV (10 44 erg/s) 




Lum 2-10 keV ( 1 44 erg/s) 



Fig. 1. (a) Centroid Fe K a line as a function of the L<2-iokeV- The radio-loud QSOs are plotted as empty triangles, 
(b) Equivalent width of the narrow line significantly detected, i.e. F-test >95%, as a function of the L2-iofceV. The 
radio-loud QSOs are also plotted as empty triangles. 



in last column of Table Q Introducing this spectral fea- 
ture significantly improved the fit statistic in 3 sources 
(0007+106, 0050+124, and 1116+215). As a further test 
we used the GAUSS model with a free to vary, instead of 
the LAOR model, and iron lines with a broad profile re- 
sulted being detected in the same objects. The intensity, 
equivalent width, and the F test for each line are shown 
in TableEl In 0007+106, an iron fluorescence emission line 
with a broad profile was detected without the additional 
presence of a narrow core, while in the last two, 0050+124 
and 1116+215, both the narrow and the broad emission 
feature are present. This limited statistic, however, does 
not allow more detailed investigation of the line parame- 
ters. 

For the QSOs with a significant Fe K a emission line we 
also tested a model with a neutral reflection componen t 
(PEXRAV in XSPEC, see iMaedziarz fe ZdziarskJ Il995|) . 
This parametrization led to no significant improvement 
in goodness of fit in the vast majority of those sources 
with resulting values of R — VI/2tt (i.e. the solid angle 
subtended by the cold reflecting material to the X-ray 
source located above it) basically unconstrained but con- 
sistent with zero. As expected, the bandpass (0.3-12 keV) 
and the limited statistics usually affecting the high energy 
(E £ 8 keV) region of the EPIC spectra are not suit- 



able to constrain the strength of this bump-like emission 
comp onent that peaked at ~ 40 keV l)Lightman & White! 
Il988h . Nonetheless, in two cases, 0007+106 and 1309+355, 
introduction of a reflection component significantly im- 
proved the fit at > 98% confidence level with a resulting R 
= 2. 61^2 and R = 2.9+-Q, respectively. Such large values 
of R indicates that simple slab geometry for the reprocess- 
ing medium is inadequate, compared to a partial-covering 
scenario with an additional spectral component emerging 
at high energies or an anisotropic geometry, which offer 
more likely explanations. The iron line parameters are not 
affected when the reflection model parametrization is ap- 
plied. 



4. Discussion 

Our study reveals that about 50% of the QSOs in the 
sample show a significant iron fluorescent emission line in 
their spectra. Iron K+shell emission therefore is a common 
feature in the X-ray spectrum of optically-selected QSOs. 
In particular, we detected 20 Fe K Q lines with a narrow 
profile and 3 broad lines. In the following section we focus 
on properties of the former feature and in Scct. l4.2l discuss 
characteristics of the broad emission lines. 
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4.1. The narrow iron line 

The energy centroid of all but four narrow lines (Fig. ^i) 
is consistent with matter in low ionization stages (Fel- 
XVII). The mean value of the energy for the 20 signif- 
icant lines detected is 6.53 ± 0.09 keV with a dispersion 
lower than 0.3 keV. When only the 16 detected low ionized 
Fe lines are considered, the value found is 6.41+q'q 3 keV 
with a dispersion lower than 0.05 keV. This best simul- 
taneous estimate of the energy of the line and its in- 
trinsic spread was calculated using t he ma ximum like- 
lihood technique of iMaccacaro et all {l988). The errors 
were obtained through the 68% confidence contour level. 
This finding turn s out t o contradict the ASCA results of 
iReeves fe Turner! l|2000|) . who found > 50% of the de- 
tected Fe lines at energies >6.4 keV in a sample of QSO 
with similar luminosities (see also lOeorge et a.1.1 2000). 
Nevertheless, considering only the common objects in both 
samples (i.e. 0050+124, 1116+215, and 1226+023) the 
only minor discrepancy is found in the case of 0050+124 
where an ionized line was only marginally detected in 
the ASCA data. The line detected in 1226+023 is not 
significant in none of the spectra, and the two centroid 
energies of 1116+215 obtained with ASCA and XMM- 
Newton are compatible within the errors. On the other 
hand, our results are compatible with other recent XMM- 
Newton studi es of QSOs; i n a sa mple of five high lumi- 
nosity RQQs. lPage et ahl l)2004ch find that the detected 
iron emission lines are mainly neutral. In a larger sample 
analysis of low-z QS Os performed with XMM-Newton, 
IPorauet et al. 1 1)2004() detect an iron line in twelve out of 
21 QSOs, and the values of the centroid of the Fe lines 
agree in all but two sources. Our finding also matches well 
the properties of iron K Q emission in Seyfert 1 galaxies ob- 
served with XMM-Newton and Chandra, where the mean 
peak of the narrow line is observed at (Er-„,) » 6.4 keV 
l)Nandra et~al"lll997t lYaaoob fe Padmanabnanll2004|) . 

The emission lines with energy E#- Q > 6.4 keV in 
1115+407 and 1116+215 are both consistent with emis- 
sion due to iron at very high ionization stages (i.e. He-and 
H-like, respectively); while in the case of the line observed 
in 1402+261, the observed energy exceeds any ionization 
states for the K Q emission. Reeves et al. (2004) have ex- 
plained the origin of this feature as a blue-shifted line 
originating in an accretion disk with high inclination. 

The intensity of the significantly detected narrow 
lines appears to be independent of the X-ray luminosity 
for L< 10 45 erg s _1 , (Fig. QJ); see the next Section for a 
more detailed discussion), with an average value (EW) 
= 80+20 eV and an intrinsic dispersion o~ew < 40 eV. 
Suc h a value is lower than 163+17, i.e. the (EW) inferred 
bv IReeves fc Turnerl |2000) for the radio-quiet QSOs in 
their sample. Nevertheless, the values for the EW of the 
de tected iron l i nes are consistent with the measurements 
of IPage et all l)2004b|) who using XMM-Newton data of 
five QSOs find EW values within the range 30-180 eV. 
The result is also compatible with the typical values found 
for Seyfert galaxies; the mean EW obtained by ASCA for 



a sample of Seyferts was 114+10 eV llNandra et alJll997h . 
and more recently, high resolution Chandra observations 
measured a weighted mean of 6 5 eV for the narrow cores 
llYanoob fe Pa.dmanabnanlEool. 

The average properties of Fe K-shell narrow emis- 
sion lines observed for our objects are therefore very 
similar to those inferred for the low-luminosity AGNs. 
Consequently, the origin of this feature in QSOs and 
Seyferts is expected to be the same. It is widely accepted 
that the fluorescence emission line emerges in 'cold' mat- 
ter distant from the central X-ray source consistent to its 
energy and narrow profile. Even if the emitting medium 
is located in the optically-thick accretion disk, regions 
near the radius of marginal stability for a not-rotating 
or spinning black hole (i.e. 6 and 1.24 gravitational radii, 
respectively) are ruled out. Furthermore, zones of the disk 
very near to the X-ray source are likely to be ionized 
given the high irradiation. A substantial fraction of the 
line flux could be due to reflection of the primary con- 
tinuu m off the inner walls of the tor us into our line of 
sight l)GhiseHini. Haardt. fe Mattil994j) . This parsec-scale 
structure is a key component of AGN Unification mod- 
els, and is indeed tho ught to subtend a large solid angle 
to the X-ray source l|Ma,tt et alj EoOO). Assuming stan- 
dard values for Type 1 AGNs (i.e. a torus column density 
Nh ~10 24 cm -2 , a viewing angle ^ 20°, and a half open- 
ing angle of the torus 9 ~ 30°), the predicted EW of the 
line produced by a torus with s olar abundances is ~ 100 
eV dOhisellim. Haardt. fc Ma.tdll99i . in agreement with 
our observations. 

Finally, we should consider the possibility of a contri- 
bution, at least partially, from the optical-to-ultraviolet 
Broad Line Region (BLR) to the line emission. On the ba- 
sis of variability and emiss ion properties (energy, width, 
EW), lYaaoob et all l|200lh inferred that the bulk of the 
iron line in NGC 5548 is indeed consistent with emission 
in the ou ter BLR. Such inte rpretation has been also pro- 
posed bv iKaspi et alJ l|2002() for the Fe K a line found in 
a very long Chandra HETGS observation of NGC 3783. 
The equivalent hydrogen column of the BLR, N#, can be 
estimated using the EW of the iron line, the index of the 
power law, and the fraction of sky covered by the clouds 
l|Yaaoob et al]l200l|) . Considering the mean index of the 
power law obtained in analysis of the QSOs in the sam- 
ple, (r 2 _ 12fcey ) = 1.87 JPiconcelli et alll2004bh. and the 
(EW)=80 eV, we get 



N 



H 



2.5 



0.35 
17 



x 10 23 cm 



(1) 



where f c is the covering factor. Assuming spherical cover- 
ing, i.e. f c = 1, the equivalent hydrogen column density 
is Nh — 9 x 10 22 cm -2 . This value can be considered 
as a lower limit to the absorption column and is, there- 
fore, in agreement with typical values of BLR clouds, i.e. 
N H ~ 10 23 cm- 2 . 

Emission lines from ionized iron (i.e. with an energy 
> 6.4 keV) can be produced by reflection in an ion- 
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ized accretion disk (e.g. iReeves et all (|200lf) ') or, alterna- 
tively, by warm scat tering material along the line of sight 
i Bianchi et aLll2004() . The latter mechanism also appears 
likely in the case of 1116+215, for which the presence of 
a warm medium is suggested by the detection of warm 
absor ber signatures in the sof t X-ray portion of the spec- 
trum ijPiconcelli et al.ll2004bJl . Similarly, the strength of 
the Fe line in 2214+139 is consistent with an origin in the 
low-ionization compone nt of the two-zone warm absorber 
observed in this source i Piconcelli et al.ll2f)04al) . 



4.2. The broad iron line 

In contrast to the narrow Fe K a emission line detected 
in 20 of the sources, the broad component seems to be 
less common in QSOs with only 3 out of 38 significant 
detections. 

Emission from an accretion disk reaching down to 
a few gravitational radii of the black hole is the 
most likely explanation for iron K a lines with a broad 
profile. Such a skewed velocity profile indeed sug- 
gests that fluorescence takes place in a strong gravity 
regime (see lFabian et al.ll2000t for a review). Except for 
a very few remarkable cases || 3 orauet fc Reeves! Hf 



IComastri. Brusa. & Oivanol 12004 ). most of the detected 
relativistic iron li nes have been det ected in nearby 
Seyfert-like AGNs l|Nandra et alJli997j) . Their presence 
in QSOs appears to be rare, which has been interpreted 
as evidence of different physical properties of the accre- 
tion disk in these objects, such as high ionization or small 
transition radius to radiatively inefficient accretion flows. 
Nevertheless a selection effect due to the poor statistical 
quality of most X-ray quasar spectra must be borne in 
mind, and detection of broad Fe lines in three quasars of 
our sample therefore appears noteworthy. 

The broad features observed in 0007+106 and 
0050+124 could also be explained as contributions from 
different iron lines at different ionization levels as sug- 
gested by the presence of some posi tive data-tct- rati o 
residuals in the range 6-7 keV (see also lGallo et al.lE'004h . 
However, current XMM-Newton data does not allow us to 
investigate this hypothesis further. Similarly, the limited 
signal-to-noise of the observation of 1116+215 hampers 
a detailed analysis of the broad iron lines and affects the 
significance of the inferred spectral parameters, i.e. leav- 
ing almost unconstrained the value of the line intensity ( 



EW=500 



-1100 



eV). 



4.2.1. Behavior of the EW with L 2 _iokeV 

The Baldwin effect was first reporte d in the o ptical and 
ultraviolet emission lines of QSOs. iBaldwirJ l(l977h ob- 
served a decrease in the equivalent width of the C/y 
emission line with increasing ultraviolet luminosity. This 
anticorrelation between the luminosity and the equiva- 
lent width was a lso found in the X-ray band for the Fe 
K Q emission line. Ilwasawa fc Taniguchil l|l993h studied 37 



AGNs observed with Ginga and found the relation EW(Fe 
K ) oc L-°f°£°y 03 . A similar study o f 39 AGNs performe d 
with ASCA als o confirmed t his effe ct l|Nandra et al.ll997h . 
More recently, IPage et alJ l|2004al) obtained a relation of 
EW(Fe K a ) oc L^-io^y 08 for a sample of 53 Type I AGNs 



observed with XMM-Newton. 

In Fig. |2 we show the relation between the 2-10 keV 
luminosity and the equivalent width of the narrow Fe K Q 
line for the objects in our sample. In order to avoid mix- 
ture with RLQs, where hard X-ray emission could be con- 
taminated by the relativistic jet components, we have con- 
sidered here only RQQs. 

Several different techniques have been compared to 
study the behavior of the EW with luminosity. 

First, linear regressions of logEW versus logL 2 _iokeV 
were performed using ASURV (Astronomy Survival 
Analysis, Fiegelson & Nelson, 1985) package. ASURV fa- 
cilitates a correct statistical analysis when censored data 
(upper limits in this case) are present, although errors are 
not considered in the calculations (see below) . This allows 
us to study the relationship between EW and L2-iofceV, 
including all significant and non-significant (the upper 
limits) detection of the neutral narrow lines. For the four 
ionized lines found in the analysis, an extra neutral compo- 
nent was added to the fit model. This component, modeled 
by a Gaussian line with er=0, is not significantly required 
in any of the four sources. The upper limits of the corre- 
sponding EWs can also be found in Tabled 

Within the ASURV package, we performed a 
Spearman-Rank test considering the detections and up- 
per limits, and found only a 30% probability that log EW 
and logL2-iokcV are correlated. The fit calculated with 
the Buckley-James method gives a power law index of 
a = -0.06 ± 0.20, where EWf F eKJ o c L%_ 1QkeV . The 
result obtained by IPage et all l|2004afl with the same 
method, a = —0.17 ± 0.08, is compatible with our value 
within the errors. 

The disadvantage of using the ASURV package in 
the data analysis is that actual detections are considered 
without taking errors into account, thereby underestimat- 
ing the uncertainties in determining the slope in certain 
cases. In order to avoid this problem, we calculated 
the mean EW in five different luminosity bins with 
a similar number of objects in each bin, therefore we 
studied the relationship between EW and L2-wkeV using 
these five bins. The weighted mean EWs were calculated 
considering asymmetric errors statistics (Barlow 2004). 
The index of the power law obtained is a = —0.07 + 0.20, 
very similar to the results given by ASURV when upper 
limits are considered in the fit. Fig. E{b) shows the data 
together with the fit. 

Fig. 13 shows a comparison of the data-to-continuum 
residuals in the 4.5-7 KeV energy band, obtained by com- 
bining the ratios corresponding to objects selected in five 
luminosity intervals. For each interval, we considered the 
best fit model of each QSO but excluded the iron line 
in cases where it is significantly detected. The unbinned 
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Fig. 2. Relation between the equivalent width of the Fe K Q line and the L2-iofceV for the radio quiet quasars. (Upper 
panel) The plot includes measurements for each set of data. The solid line shows the fit to all the detections while 
for the dashed line only the significant detections are considered. (Lower panel) The plot represents all measurements 
binned in five groups. The solid lines show the fit to the five points. 



data-to-model ratio of each QSO is translated into the 
rest frame and grouped in bins of 0.2 keV size. Each plot 
also includes a Gaussian emission line broadened to the 
intrinsic resolution of the pn at the energy of the iron line, 
i.e. FWHM=150 eV. The plots show residuals compatible 
with the neutral iron line being narrow. The residuals in 
the last bin are larger and randomly distributed. 

An anticorrelation between EW and hard X-ray 
luminosity is also compatible with our data. However, 
other behaviors different from a power law cannot be 
ruled out. The data actually suggest a constant EW 
for L 2 -io kov ~ 10 45 erg s _1 and a stepwise decrease 
above this threshold. As can be seen in Fig. E{b) the 
anticorrelation between EW and L2-iokoV is actually 
determined mainly by the highest luminosity bin. 

The physical origin of variation in the iron line 
equivalent width with luminosity is related to where 
the line is originated. Assuming that the Fe K Q line is 
produced by fluorescence, vario u s scen arios have been 
proposed. Ilwasawa fe Tanieuchil (|l093^) suggested that 
the emission line is produced in the clouds of the Broad 
Line Region. The X-ray primary emission ionizes these 
clouds explaining the ob served decrease in the equivalent 
width with luminosity. iNandra et al.l l|l997h proposed 
that for iron lines, which are most likely produced in 
the accretion disk, an increasing degree of ionization 
with luminosity could also account for the Baldwin 
effect. Assuming the em ission line is originated in the 
torus. IPaee et alJ l|2004a|) suggested that increasing X-ray 
emission could cause a decrease in the covering factor of 
the torus and hence a decrease in the equivalent width 
of the iron line. A step function behavior of the value 
of EW with luminosity favors our explanation in terms 
of the increase in the degre e of ionization suggested by 
Ilwasawa fc Taniguchil l|l99. c ifl . When the material reaches 
a particular ionization state, the emission line could be 



suppressed by resonant trapping and/or full electron 
stripping and the EW would then abruptly decrease. In 
such a case, the material where the line is emitted should 
be most likely identify with the BLR or accretion disk. In 
the last case, the average unresolved profile indicates an 
origin in the outermost regions of the accretion disk. 



Nevertheless, our analysis cannot rule out a smooth 
power law dependence , EW(F eK rt ) cx L%_ lokeV , as sug- 
gested by IPage et alJ l|2004a|) . The slope of the power 
law ranges between -0.27 and 0.14, taking into account 
the different analysis and their corresponding uncertain- 
ties. More highly luminous QSO with high signal-to- 
noise need to be added to this analysis to reduce large 
uncertainties in determining the slope. A possible diver- 
gence from previous works could be due to a selection 
effect, as they claim the presence of X-ray Baldwin ef- 
fect and include both RQQs and RLQs in their sample. 
In particular, most of the objects with high luminosity, 
Ly-inke.v > 10 45 erg s -1 , were r adio loud. However, as 
shown bv lReeves fe Turner! l)2000|) . there is a strong anti- 
correlation between iron K a EW and radio loudness. This 
hypothesis is also strongly supported by anticorrelation 
between radio loudness an d the strength of the s oft ex- 
cess component studied in IPiconcelli et all l|2004bf) . This 
trend could be due to the likely presence of an additional 
jet emission component that weakens the strength of the 
iron line. Inclusion of RLQs in the highest 2-10 keV lumi- 
nosity bins in previous works could therefore have affected 
the results on the X-ray Baldwin effect. The choice of our 
sample minimizes this bias. 

A similar analysis of the radio loud objects is ham- 
pered due to the low number of RLQs in the sample. 
Furthermore, a significant Fe K Q line was detected in only 
one out of five RLQs. 
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Fig. 3. Mean data-to-model ratios in five luminosity bins, assuming for each QSO the best fit model and excluding 
the iron line and correcting by redshift. The observed lines represent a Gaussian profiles with centroid energy 6.4 keV, 
and width equal to the intrinsic pn resolution, normalized to the 6.4 keV data point. 



5. Summary and Conclusions 

We present here an analysis of the K a fluorescent 
iron emission line in the XMM-Newton spectra of 38 
PG QSOs, including the presence and properties of 
emission lines with both narrow and broad profiles. 

A narrow unresolved (i.e. FWHM< 7000 km s _1 ) 
emission Fe K Q line was detected in ~ 50% of the 
objects studied. The energy centroids indicate that the 



majority of the lines originated in cold material (Fe 
I-XVII). The mean equivalent width for the narrow line 



is (EW) 



SO 



+30 
-20 



eV with a dispersion oew < 40 eV. 



Such properties are consistent with several scenarios for 
the origin of the line. Fluorescence likely occurs in the 
outermost parts of the accretion disk and/or in even 
more external regions, such as the molecular torus or the 
clouds of the Broad Line Region. For two objects, i. e. 
1115+407 and 1116+215, the line energies are associated 
with the highest ionization stages of iron, suggesting 
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emission from an ionized accretion disk, due to scattering 
by warm material located in the line of sight or due to 
the presence of an outflow. 

We also investigated the presence of an X-ray Baldwin 
effect for the Fe K a line of the RQQs. The RLQs were 
excluded in order to avoid possible contamination from 
a relativistic jet component in the X-ray emission, 
which could have biased the anticorrelations published 
in the literature prior to our study. However, no clear 
correlation was found between strength of the iron line 
and 2-10 keV luminosity. Assuming a power law behavior 
EW(FeK a ) oc L2— 10 fee V' a turned out to in the range 
between -0.3 and 0.14. 

Broad relativistic lines were detected in three objects, 
less than 10% of the QSOs in the sample. These objects 
have low luminosities (L 2 -iokeV ~ 10 44 erg s _1 ) suggest- 
ing a possible anticorrelation between the X-ray lumi- 
nosity (hence the ionization state of the uppermost disk 
layers) and the strength of lines originating in the disk, as 
has b een proposed by iNavakshin. Kazanas. & Kallmanl 
(200(J. However, a selection effect due to limited statistics 
affecting the most distant (very luminous) QSOs prevents 
us from drawing firm conclusions on this point. 

After comparing our results with those of recent 
XMM-Newton and Chandra observations of Seyfert 
galaxies, it seems that characteristics of the Fe K-shell 
emission lines detected in the spectra do not significantly 
change as a function of the X-ray luminosity. We report 
average values of the EW and the energy centroid of the 
line which are both consistent with tho se in Seyfert galax- 
ies. Moreover, IPiconcelli et alJ l|2004bJ ) also found a sim- 
ilarity between the properties of the cold and warm ab- 
sorption in QSOs and Seyfert galaxies. Our analysis there- 
fore allows an extension of the well studied observational 
properties of the nearby X-ray bright Seyferts to the QSO 
realm, at least for those optically selected. 
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